The current phase relation in Josephson tunnel junctions. 
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The J(ip) relation in SFIFS, SNINS and SIS tunnel junctions is studied. The method for analytical 
solution of linearized Usadel equations has been developed and applied to these structures. It is 
shown that the Josephson current across the structure has the sum of sin ip and sin 2tp components. 
Two different physical mechanisms are responsible for the sign of sin2yj. The first one is the 
depairing by current which contributes positively to the sin2y> term, while the second one is the 
finite transparency of SF or SN interfaces which provides the negative contribution. In SFIFS 
junctions, where the first harmonic vanishes at "0" - "n" transition, the calculated second harmonic 
fully determines the J(<p) curve. 



PACS numbers: 74.50.+r, 74.80.Dm, 75.30.Et 



It is well known that tunnel SIS Josephson junctions 
have sinusoidal current-phase relation, while with the de- 
crease of the barrier transparency deviations from sin <p 
take place (seei^ for the review). The sign of second har- 
monic is important for many applications, in particular 
injunctions with a more complex structure like SNINS or 
SFIFS, where N is a normal metal and F is a weak metal- 
lic ferromagne1p2i2i4. To analyze this problem selfconsis- 
tently, one should go beyond the approximation which is 
usually used and is called "Rigid boundary conditions" 
(BBC) . 

The RBC method is an effective tool extensively used 
earlier for theoretical study of the proximity and Joseph- 
son effects^. This method is based on the assumption 
that all nonlinear and nonequlibrium effects in a Joseph- 
son structure are located in a "weak link" connecting 
two superconducting electrodes. The back influence of 
these effects on superconductivity in the electrodes is ne- 
glected. The RBC are valid if a junction has the constric- 
tion geometry. The quantitative criteria for the validity 
of RBC for planar SIS tunnel junctions, SS'S sandwiches 
and variable thickness bridges were studied only numer- 
ically for some parameter ranges^. The main technical 
difficulty in formulating the analytical criteria of RBC 
validity is to find the solution of equations describing the 
perturbation of superconducting state in S electrodes. In 
this paper we will attack this problem by finding the solu- 
tion of linearized Usadel equations^. We will also use this 
solution to formulate the corrections to previous results 
obtained in RBC approximation. 

The junction model. Let us consider the structure 
of SFIFS type, where for simplicity the parameters of the 
SF bilayers are equal to each other. We assume that the 
S layers are bulk and that the dirty limit conditions are 
fulfilled in the S and F metals. We assume further that 
F metals are weak monodomain ferromagnets with zero 
electron-phonon interaction constant and the FS inter- 
faces are not magnetically active. We will restrict our- 
selves to the case of parallel orientation of the exchange 



fields H in the ferromagnets. The results obtained for 
SFIFS junctions cross over to SNINS and SIS in corre- 
sponding limits. 

Under the above assumptions the problem is re- 
duced to the solution of the one-dimensional Usadel 
equations^ in S- and F-layers and matching these so- 
lutions by the appropriate boundary conditions^ We 
choose the x axis perpendicular to the plane of the inter- 
faces with the origin at the central barrier I and introduce 
indexes L (left), R (right) and I for description the ma- 
terials and interfaces parameters of the SFIFS structure 
located on the left and right sides from the central barrier 
and at this central barrier, respectively. 

The Usadel functions G and F obey the normalization 
condition + F^Fl^ = 1, which allows the following 
paramctrization in terms of the new function 



G^j — 



<J> <J>* 



F„, 



(f> <f>* 



(1) 



The quantity uj = uj + iH corresponds to the general case 
when the exchange field H is present. However, in the S 
layers H = and we have simply u) = u. 

The Usadel equations 5 in the S and F layers have the 
form 



nT c d 
ujGs dx 
^ 2 ttT c d 
F ujGf dx 



$ s = 

- $p 



-A, 
0, 



(2) 
(3) 



where G w = uj/^/uj 2 + <f> u <f>*_ u , 2 = w + iH in a fer- 
romagnet (H is the exchange field), ui — to in S and 
N metals, T c and A are the critical temperature and 
the pair potential in a superconductor, lo = irT(2n + 1) 
are the Matsubara frequencies and £s(F) are the coher- 
ence lengths related to the diffusion constants Ds(f) a s 
£s(F) — \/Ds(f)/^ 7T T c . The pair potential satisfies the 
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self-consistency equations 



and at x = ±cIf 



Aln- +7 rT £ - = 0. (4) 

UJ— — OC 



In the case of SFIFS tunnel junction in quasi-one di- 
mensional geometry the boundary conditions at the junc- 
tion plane (x = 0) read 
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G F,L d 



1BI- 



UlL dx 
£fGfl,r d 



F,L 



FL,R 



cjfl ax 



±G 



<9x \ wr 

with 7s/ = R N Ai/p F ^F, 



F,R, 
$F,L 



(5) 
. (6) 



where the indices L and i? refer to the left- and right- 
hand side of the junction, respectively, i?jv and Ai are 
the normal resistance and the area of FIF interface. 

The boundary conditions at the SF interfaces (x = 
=FcZf) have the formi 



tsG 2 s tk d x £,FG 2 Fk d 



dx 



S.k 



4> 



F,k, 



LJ k dx 

<f>S,k <f>F,k 



, £,FG F ,k d 

±7s — TT®F,k = G S ,k 

U! k OX \ U> Ul k 

with 7s = RbA b /pf£,f, 7 = Ps£s I 'pf£f, 



(7) 
(8) 



where Rb and Ab are the resistance and the area of the 
SF interfaces; ps(F) is the resistivity of the S (F) layer; 
k = L. R.. Both of these conditions ensure continuity of 
the supercurrent. 

We will also suppose that due to low transparency of 
the FIF interface the Josephson current is much smaller 
that the depairing current of superconducting electrodes 
so that the suppression of superconductivity in the inte- 
rior of the electrodes can be neglected and at x — > ±oo 



1$ 



S,k 



= A , 



(9) 



where Ao is the magnitude of bulk order parameter. 
The limit of small F layer thickness. In this limit 



Ak = A ,k =f 7b 



G 



F.k 



Gs,k + Uk^BM /7TT ( 



■Bk 



(13) 



.4, 



^>R,L^S.kGs,i 



dp . , 

o,fe = —77; — T^FT' 7bm = 7b t - - (14) 

<*> {Gs,k + uj k 7BM/^T c ) £p 



Expression valid if 75 -C 7b/. Substitution of ltTT|) 
and 113|l into the boundary condition at x = leads 
to 



d Gp k w ^Gp u 



'u k G 2 



(15) 



S.k 



where 7a/ = jdp/^p and reduce boundary problem 
© to the solution of equations J2J, Q in the S-layers 
with the boundary conditions ©, l|15l) . At H = and 
(lBid/£,F) ^> 1 expression ( I15|) reduces to the known 
result for SN bilayer^ 

The linearized Usadel equations. Following RBC 
approximation we will start with the assumption that the 
suppression of superconductivity in S layer is weak and 
the solution of Usadel equations in the superconductor 
has the form 



Gs,k 

G\h 



A , fc 

Go + G\.k 
GV, 



$ 1)fe , A = A , fc + A ljfe , 



Go 



A 2 
^0 



(16) 
(17) 



A5 ifc $ 1|fc + Ao lfc $I ifc 



c^ 2 + A 2 



where Ao,fc = Ao exp {±iip/2 + illx/^s} , f is order pa- 
rameter phase difference across the barrier and the coef- 
ficient U describes the linear growth of phase difference 
due to the supercurrent in the electrodes. Corrections to 
Ao and &s,k are supposed to be small 



dp <C min £ F 




(10) 



the gradients in (|3J) are small and in the second approxi- 
mation on df/C-F the solution of 10 has the form 



-r . . D x x 2 uj k A k 
®F,k = A k + B k — + — 2 , 

4F 2 7T1 c ^ pLrF,k 



r 2 - 



(11) 



Integration constants A and B in (|ll|l can be found from 
boundary conditions at x = 



Gf,L t-j G 2 FR Gf,lGf,R 
— -£>L = — &R — 



ML 



UR 



IB I 



i^-^) (12) 

UJR LU L 



|A lifc |<A , |$i, fc |<A . 



(18) 



The approximation is valid if the right hand side of 
Ea.pSjl is also small, so that 



i— 1 / \ wGF0,kA() k 

r, k {uj) = ±7m — 1777-^2 h 7" 



'F0,fc 



ttT c G 2 



G 



ivd k 



FO.k 



V^ + A 2 G 2 ' 



(19) 



(20) 



where = (G + uik^BM /ttTc) > and |S(w)| <C A . From 
the structure of the linearized Usadel equations and the 
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boundary conditions i|19[l it follows that there are first 
order corrections only to the magnitudes O and Ai of 
functions 4>i and A^/. respectively, while the phases of 
all of these functions coincide with those of Ao,&. In this 
case 

$i, fe = 9exp{±i|}, Ai, fe = Aiexp{±i|} (21) 

and due to the symmetry of the structure we have 
lor = u>l = 5, Gfqm = Gfo, "Ok = ■& 



A),fc n /, .<P\ n 
_ = 0,exp|±z-}, Co 



uG 



The solution of flfy . is 



OO j 

Ai = 2^ °n cxp(-q — — — ), 



is 



(31) 



Sny/uj 2 + A 2 , x-d F 

9 = / o a =t — ^r ex P(~9o— ^— ), 



where the coefficients and satisfy the equation 

q 2 n nT c 



S 2 =7TT^ 



, (32) 



Gfo 
~Gn~d 



to w 

±1M *E C0S 2 + 



+ l [ lM—^r + 2 Gfo sin — 



(22) 
(23) 



To write Q22JI , we also used the fact that in the first order 
with respect to |H(w)| the magnitudes of functions &s,k 
in 1|13|) equal to Ao and that Gs = Go- 

Substituting lfTS)l. lf?T|l into @, ©, we arrive at the 
following boundary problem for and Ai 



-e s 



7TT n 



d 2 



y^ 2 + A 2 dx 2 



e + e = Ai, 



(24) 



E 



= A P(^) 

' (x/^+A 2 - 7rT c g2 ) V^+A 2 



(33) 



and P(ip,uj) = ReSfl(w) cos(</?/2) + ImHi{(£j) sin(<p/2). 
Multiplying Eq.JSJl on w 2 (u> 2 + A 2 ,) -3 / 2 , summing both 
sides of this equation on lj and making use of ll-il'l) one 
can transform ll.'j.'ill into the system of equations for the 
coefficients Sq which yield 



TT^AottV 



(34) 



where 



Ai 



In 



T 



oo 1 
^ \f,i\ 



0J — — OO 



LJ — — OO 



(25) 



A(n,^) = 



jmKAQ) + -J-K 2 (n)(l - cos^) 

1BI 



d 

£s^0(±d F ) = ReS fe (w)cos- ±ImS fe (w)sin- 



(26) 



6(±co) = 0. 



(27) 



Due to the symmetry of the problem it is enough to 
solve the equations l|24|) - (|27|l only in one of the electrodes, 
namely, for x > (If- Using the equation for Aq(T) 



Iu^+ttT J- 

T r ^ \Uj\ 



1 



UJ — — OO 



(28) 



and the symmetry relation Q(u>) = Q(—lo) we can rewrite 
the selfconsistency equation in the form 



^ t:T c lj 2 2 8 2 
= n ^ (lu 2 + Aq) 2 'dx 2 



UJ>0 



A 2 



^ (w 2 + A 2 ) 3 / 2 ' 



(29) 



(30) 



pG + (Hq+pty-fBAi/nTc 
G (p 2 + q 2 ) 



H f2 

q = 2^ BM — (j BM — + G ), (36) 



_ jj2 

p=l+ , „. 2 71m + 2G ^r7BA/- (37) 



f2 

' 7rT r 



Here SI = ixT(2m + 1) are the Matsubara frequencies. 

As a result, the solution of the boundary problem 124(1 - 
(|2*7Jl has the form 



Al = _ E S2( » 2 + a 2 ) 2 A ^ 

n>o y 0J 

(38) 



4 



= -nTj2 



n>o 



TrT c A n 2 q n A(n, if) exp(-go^) 
E 2 (fi2 + Ag)2(l-^r c g2 Go / w ) ' 



(39) 



In particular, at x = d F from (|38[1 and l|39|) we have 

= - 7 mSfx - — S F2 (1 - cos^), (40) 
A 7b/ 



and Eq. (|46|l we can write down the supercurrent / across 
the SFIFS junction in the form 



where 
Jo 



/ == (Jo + Ju) simp + J12 sin2<^, (48) 



1 t m 



eR N ^ u 2 + C 2 A 2 



j0 



n>o 



^ 2 + A 2 ) 2 (l ~ ttT^Go/uj 



(41) 



Ju = - 



2ttT 

cRn 



E 



UJ — — QO 



^O^O 



(£ 2 



G 2 A 2 ) 2 



7m — t Zjfi~ 



■>F1 



7rT c 2 g n ^ 2 (0) 



o>o 



S 2 (fi 2 + A 2 ) 2 (l-^T c(Z 2 Go/c)- 



(42) 



To calculate the sums (|41[1 , ()42[) one needs to know the 
expression for the coefficients go which can be in gen- 
eral obtained from numerical solution of Ea. i|32|) . Since 
the main contribution to the sums (|41|l , (|42(l comes from 
large fi, the asymptotic behavior of qn at large f2 can be 
used 



2 

In = " 



7rT R 



(43) 



The developed method is valid if the following condi- 
tion is fulfilled 



7 - 
7b/' 



! max < 1 , In 



H 2 



min {7lM.7lf}( 7rT ) 



«1, 
(44) 



7b < 7b/- 

Therefore for the function &s,k in Ea. (|14|l we get 

$s,fe = (A + 9{d F )) exp , (45) 

and substituting Q45JI into Ijl3|l we finally obtain 



Ai 



7b wG GfoAo . <p 

T2« ^ sin-, 

7b/ u>V 2 2 



Go exp {±iip/2} =f (46) 
(47) 



/i = 1 + GouJ-fBAi/nTc 

From the structure of coefficients ^4_r,l we see that the 
corrections to the supercurrent across the SFIFS tunnel 
junction leads not only to the reduction of the critical 
current of the structure, but also to changes in the 3 s (ip) 
relation. 

The Jg(y) relation. Using the the standard expres- 
sion for the supercurrenljii, the boundary condition © 



7b u 2 Gfo . 7 ojcoCofi 



Ibi 



Ju 



eR N 



OO 

E 

■J — — OO 



7b/ 



A 2 r 3 

^0^0 



An 



(50) 



(£ 2 + G 2 A 2 ) 2 



7 B GfqAqGq 

7b/ i? 



7 WW/i.Si?2 



(52) 



7b/ A 



Expression (|49|1 has been obtained previously i n 9 11 . 
The ip— independent correction to it, Ju, is negative and 
describes the suppression of sin ip component of the su- 
percurrent. The first term in Eq. (|50|l proportional to 
7m takes into account the suppression of superconduc- 
tivity in S electrodes due to proximity with thin F layer. 
The last two terms proportional to 7^) describe the sup- 
pression of superconductivity by the current across the 
junction. The larger 7b and 7 the weaker is the super- 
conductivity induced into F layer and the stronger is the 
influence of this effect. 

The sign of the second harmonic Ji 2 depends on the 
relation between 73 and 7. At 7b = it is positive and 
J(<p) relation l|48|) has a maximum at ip = (p max < ir/2 . 
Such a shift was predicted earlier near T c for SIS tunnel 
junctions and is due to the suppression of superconduc- 
tivity near the barrier by a supercurrent^. Increase of 
7b leads to additional phase shifts at both SF interfaces 
and provides the mechanism for the shift of the y ma x 
into the region <p > ir/2. As a result, at sufficiently large 
7b the amplitude Ji 2 changes its sign and <^ max shifts to 
ip > 7r/2. Such a competition between suppression by a 
supercurrent and by proximity effect was first analyzed 
in the SNS junctions^ at T w T c .This fact is in the 
full agreement with the results of numerical calculations 
summarized in 2 . 

The physical reason for different signs of Ji 2 can be 
easily understood if we consider the two cases separately. 
Suppose first that 7b is finite. In this case the SFIFS 
structure may be considered as a system of three Joseph- 
son junctions in series as shown schematically in Fig.l. 
For rough estimates one can assume that the phase x of 
&F,k does not depend on w. Demanding the equality of 
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FIG. 1: The phase distribution in a SFIFS junction 



FIG. 2: Depairing by current near the tunnel barrier 



the currents across FIF and FS interfaces and taking into 
account that Iq oc 7^) <C Ici ■yg 1 for \ we will have 



X = W2 



sin 2%. 



Substituting this x m to the expression for the supercur- 
rent across FIF interface, we get 

I = Ic sin(< ( 5 — -2- sin ip) « Ic(sm.(f — sin 293). (53) 
Ici Ibi 

Therefore with increasing 73 the phase partly jumps at 
the FS interfaces leading to a continuous crossover from 
the Josephson effect lumped at x — to the phase drop 
distributed at \x\ < dp. In a full agreement with the 
theory of double barrier devices^ this crossover results in 
appearance of second harmonic in Jsif) with negative 
sign which provides maximum Js(ip) achieved at ip > 
tt/2. 

If 1b = 0, the structure is always lumped at x = and 
the main effect is the suppression of superconductivity 
by supercurrent in the vicinity of FIF interface as shown 
schematically in Fig. 2. The resulting contribution to the 
full current is 



2 w 



■) sin<p. 



T 1 / a y 99, . A _ , 

4, oc 7b/(Ao - £s — )sm</? ex (1 : 

OX JBI IBI 

(54) 

It follows directly from l|54|l that the amplitude of the 
second harmonic is positive. 

The competition of the above two mechanisms of I(tp) 
deformation is clearly seen from Eq. (|51|l . 

The general expressions (|49H - <|51|) can be simplified in 
several limiting cases. 

In the symmetric SNINS tunnel junctions H = in 
both electrodes and in the first approximation from i|49|) 
the earlier result fromi^ is reproduced 



Jo 



2ttT 



eR N ^ 



E 



A 2 
^0 



A§)6H 



(55) 



Q(u) = (1 + 2G u 1b m/kT c + (lu 1bm /ttT c ) 2 ) 



while 1)50(1 and H51[) reduce to 
AttT 



J 11 



eRN 



7A/S4 + — — £5 H — E6 

IBI JBI 



Jv 



2ttT 
cRn 



1b v 7 v 

2^7 Ijq 



IBI 



IBI 



where 



E 

oj>0 



E 



(^ 2 + A2)e 2 M' 



A 2 tf 2 



LU>0 



(w 2 + A2)9 5 / 2 H' 



(56) 
(57) 

(58) 
(59) 
(60) 
(61) 

jj- t l±q)-<o~> -\ui) 

and Go = u>/ \Joj 2 + A 2 ,. 

In the limit 7 — ► 1, ff,7M>7B)7sM — * the SFIFS 
structure transforms into SIS tunnel junction. In this 
case 

Co = 1, = I A o + 9(<*f)] cxp {±*^/2} , (62) 

7rT c A f> 2 g o sin 2 f 



E 

w>0 



goAg^S£2 

(^ 2 + A2)e 2 H' 



E 



A 4 

' (c^ + a 2 )^ 5 / 2 ^)' 



e(d 



;^E 



75/ s ^ E 2 (fi 2 + A 2 ) 2 (l- 7 rr c g 2 GoM' 

and for the supercurrent I in the first approximation 
we have the well known result of Ambegakaokar-Baratoff 
theorjii^ 



I 



2ttT ^ A 



E— — sm tp. 
,,,,, _w 2 + A§ y 



(64) 



w>0 
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Using ilML'l for Jn and J 12 it is easy to get 



Jn 



<^3, ^12 — —p. — ^3, 



eR 



N 



4 x - A ir 

3 = ^ T t WTW^' 



(65) 



(66) 



and the full current across the tunnel junctions is 
r A 



- tanh — - 2E 3 
2 2T 



, A £ 3 . , . 

smyH — — sin Zip. (67) 



The critical current achieves at phase difference ip c 



7T 4£ 3 a A 

¥V = tanh — , 

Yc 2 7T 2T' 



(68) 



and equals to 



A 



- tanh — - 2£ 3 
2 2T 13 



and at T — > the I(ip) simplifies to 



A 7T 1.92 /ttT ( 



ei?jv 2 7b/ V Ao 



3/2' 



(69) 



At T w T c Eqs. lfH5fl transform to the result obtained ir*i&. 

Conclusions. In summary, we have studied the 
current-phase relations Js(<p) in SFIFS, SNINS and SIS 
junctions in the regime when the second harmonic of 
Js(f) is not small. To solve this problem selfconsistently, 



we have developed the analytical method for solving the 
linearized Usadel equations. This solution describes a 
weak suppression of superconducting state in a supercon- 
ductor caused either by proximity with normal or ferro- 
magnetic material or by a current in composite SN or SF 
proximity systems. The method is rather general and can 
be applied to a wide spectrum of proximity problems. 

We have demonstrate that the full current across the 
structure (|48|l consists of the sum sin ip and sin 2ip com- 
ponents and have calculated the amplitudes (Jo + -hi) 
and J12 of these components. In SIS and SNINS struc- 
tures the corrections Jn and J12 to the previously cal- 
culated critical current Jo are small. The J((p) curve is 
slightly deformed so that the maximum value of the su- 
percurrent achieved at phase difference <p c which can be 
smaller or larger ir/2 for positive and negative sign of J12 
respectively. In SFIFS junctions Jo = at the point of 
the transition from "0" to "7r" state. It means that in 
this case the calculated values Jn and J12 determine the 
J(<p) curve. Since the amplitudes Jn and J12 may have 
comparable magnitude, the J(ip) measured experimen- 
tally can be essentially different from sin ip. The validity 
of the developed approach is determined by inequalities 
(I44|l and 75 < 75;. These conditions also determine the 
validity of rigid boundary conditions in the models^ de- 
scribing the properties of SFIFS, SNINS and SIS tunnel 
junctions. 
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